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A microscopic mechanism of low-pressure insertion and separation of chain-like 
molecules in host matrices is proposed. It is shown that the intramolecular correlations 
combined to appropriate host activities are responsible for a low-pressure condensation 
of chain molecules. This allows recover a fine structure of the isotherms and to explain 
recent experiments on the insertion of C2H2 and CO2 guest species. We argue that 
the mechanism should be dominant in low-dimensional host geometries, where the en- 
tropic effects are strongly suppressed and the major factors are the chain connectivity 
and packing. 



I. INTRODUCTION 



Insertion of guest species into host matrices 
has recently attracted a considerable interest as 
a designing tool for hydrogen storage systems [1], 
rechargeable batteries and superconducting de- 
vices [2] . Recent progress in synthesis of new com- 
posite materials [3,4] opens novel perspectives in 
separation and storage of larger molecules (e.g., 
acetylene, methane or even longer hydrocarbons) 
which have important technological applications. 
In particular, a highly controlled acetylene accom- 
modation in a metal-organic microporous material 
has recently been reported [5]. It has been found 
that due to formation of regular hydrogen bonds 
the guest molecules are arranged into dense one- 
dimensional strips at very low relative pressures. 
This allows for a safe storage at a density much 
higher than the compression limit of free acety- 
lene. In contrast, CO2 isotherms exhibit a much 
weaker adsorption in this pressure domain. This 
offers an efficient way of separation of these two 
compounds. In this context it has been empha- 
sized [6] that the fabrication of high-capacity ad- 
sorbents, sensitive to particular guest species, re- 
quires a better knowledge of the host-guest and 
guest-guest interactions. In particular, our under- 
standing of a microscopic mechanism, that allows 
to reach high densities inside a host marix at very 
low bulk pressures, would be quite helpful for a 
designing of novel storage and separation systems. 

In this Letter we report on a possible micro- 
scopic mechanism, responsible for a remarkable en- 
hancement in the insertion of chain-like molecules 
into crystalline host matrices at low bulk pressures. 
The main idea is to note that in the host-guest con- 
tact the guest acts as a thermodynamic subsys- 
tem that has an impact (energetic and entropic) 
from its bulk. This induces a correlation for the 
guest species, such that their accommodation is 
not driven only by the accessibility of the matrix 
space but also by effective (density-dependent) in- 
teractions, arising from the packing and the chain 
connectivity effects. For planar geometries this 



mechanism has already found experimental evi- 
dences [7,8] and interesting kinetic interpretations 
[9]. 



II. INSERTION MODEL 

The host is modeled as a lattice of adsorbing 
sites with the lattice spacing d and coordination 
number q. The lattice is assumed to be absolutely 
rigid, such that possible guest-induced distortions 
[10] or restructurings [11] are excluded. The guest 
molecule is modeled as a chain (length m) of hard- 
sphere segments [12] (diameter a). In this way we 
take into account the chain connectivity and the 
excluded volume effects. Therefore, the guest fluid 
is characterized by its bulk density p (or packing 
r\ = 7rpCT 3 /6) and the chain length m. The latter is 
the number of hard sphere segments which mimic 
the functional groups (e.g. m — 2 for C2H2). The 
bulk equation of state can be approximated as [12] 
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where Phs is the hard sphere contribution and 
/3 = l/(kT) is the inverse temperature. In what 
follows this equation of state will be used to make 
a link between the bulk density p and the pressure 
P. 

The adsorption potential U s (ri) is given by the 
sticky site model [13,14] 

exp(-/3£/ s (r 4 )) = 1 + A^<5(r 4 - R m ) (2) 

R m 

where is a chain segment position, R m is a po- 
sition of a lattice site and A is the stickiness pa- 
rameter. The latter can be related to the segment 
adsorption energy e through A = exp(/3e). This 
singular one-body potential allows one to perform 
an exact integration in the partition function 
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where Z re f is the partition function for the same 
system but without the adsorbing potential, and 
p r ^' (Ri, R„) is the n-body distribution func- 
tion taken at positions of the lattice sites. In gen- 
eral, even in the absence of a specific adsorption, a 
periodic host structure induces a density redistri- 
bution and particle correlations around the lattice 
strands. Nevertheless, in a low-pressure regime (on 
which focus), these effects can be neglected and the 
reference state can be approximated by the bulk 
state. Therefore, we have in eq. (3) an infinite se- 
ries on A including the correlations of all orders for 
the reference state. If only pair correlations are 
important then the problem can be mapped onto 
the lattice gas model [13,14] 

S = ZjZ reS = ]T cxp(-/?H LG ) (4) 
{U} 

with the Hamiltonian 

Hlg = W{K ^ R i)***i _ J2 (5) 

ij i 

in which {U} is a set of occupation numbers. The 
effective chemical potential /Li(Rj) and pair inter- 
action W(Ri,Hj) are closely connected with the 
properties of the fluid in the bulk phase. Namely 
[13,14], 

/? M (R;) = ln(A^ e/ (Ri,<7/2)) (6) 



Here the first term corresponds to the hard sphere 
packing effects which are dominating with increas- 
ing r\. The second term describes a depletion with 
increasing to due to the chain-chain repulsion - 
the intermolecular correlation. The last term in 
(72(c) is responsible for the intramolecular correla- 
tion inside a chain. The latter is dominating at low 
densities when the chain-chain contacts are rather 
rare and the segments correlate mainly through the 
chain connectivity. That is, a probability to find 
two randomly chosen segments in the same chain 
is rather high. The singularity at 77 — > is a flaw of 
the Percus-Yevick approximation, which, however, 
does not affect our conclusions. The coefficients 
K and Q control the corresponding magnitudes, 
which depend on the chemical nature of molecules 
(e.g. chain flexibility) and on the dimensionality of 
space (in ID, for instance, the intermolecular cor- 
relations are realized only through the chain ends). 
Note that the chain length to is chosen as an ex- 
ternal variable. This means that we do not con- 
sider the clustering process itself, assuming that 
the guest cluster composition remains fixed in the 
course of insertion. This assumption is valid if the 
adsorption is not dissociative or if the adsorption 
rate dominates the rate of dissociation in the fluid 
bulk. 



A. Effective interaction 
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where p[ e (Ri) and g™ (Rj,Rj) are respectively 
the one-body and pair correlation functions for the 
reference state (i.e., a state without specific ad- 
sorption). The latter is independent of the ad- 
sorption site position because the lattice is trans- 
lationally invariant. If a ~ d then we deal with the 
contact values of the strand-segment density pro- 
file pi(a/2) and the segment- segment correlation 
function .92(c). The contact density is approxi- 
mated by its hard wall-chain value [14] 
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Here the second term describes a chain depletion 
compared to a purely monomeric case (to =1). At 
low pressures (bulk densities) this effect is marginal 
and pi(a/2) can be replaced by pa 3 /m in this do- 
main. In the Percus-Yevick approximation the pair 
correlation is given by [15] 
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As we have already noticed, due to correlations 
in the reference state (bulk), the adsorbed guest 
species experience an effective pairwise lateral in- 
teraction /3W = — In 32(0"), where 52(c) ~ 92(d) 
is given by eq. (9). In Figure 1 this quantity is 
analyzed as function of the packing 77 at differ- 
ent chain lengths to. It is seen that (3W is a 
non-monotonic function. In a low-density (pres- 
sure) domain, dominated by the intramolecular 
correlations, the interaction is strongly attractive. 
This effect increases with the chain length. At in- 
termediate densities the chain-chain competition 
for the available space (the intermolecular corre- 
lations) leads to an effective repulsion at neigh- 
boring lattice sites. Finally, at high densities the 
chain connectivity plays only a marginal role since 
a closely packed chain configuration becomes iden- 
tical to that of disconnected hard-sphere segments. 
This leads to an effective attraction [13,14], simi- 
lar to the depletion forces in colloidal systems. Al- 
though our analysis is based on the approximate 
result eq. (9), our conclusions should hold in gen- 
eral because, from the experimental point of view 
[7-9], the density/concentration regimes (e.g. di- 
luted, semidiluted and plateau regimes) are clearly 
distinguished. 



B. Insertion isotherms 



Having determined the inputs to the lattice gas 
Hamiltonian we can find the mean-field approxi- 
mation for the free energy, which gives the adsorp- 
tion isotherm [13,14] 
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where 9 — §2nti)/N is the coverage - the guest 
density inside the matrix. Equation (10) has been 
solved numerically in order to determine the cov- 
erage 9 as a function of density r\ and chain length 
size m. Then a pressure dependence can be ob- 
tained through the equation of state (1). The re- 
sult for a dimcric fluid (m = 2) on a four-fold coor- 
dinated lattice (q = 4) is plotted in Figure 2. It is 
seen that the isotherm shape is coherent with the 
density (pressure) modulation of the effective in- 
teraction PW (see Figure 1). Namely, at low pres- 
sures when the chain-chain contacts are negligible, 
the adsorption of one bead promotes the adsorp- 
tion of the entire chain through the intramolecular 
correlations. As a result we have a remarkable in- 
crease of the guest density 9 in a very narrow pres- 
sure domain. This cooperative effect has been con- 
firmed experimentally [7]. At intermediate pres- 
sures the chains compete for the adsorption sites 
(the intermolecular correlations) and the isotherm 
exhibits a "slower" increase with a characteristic 
inflection point. The latter signals a cross-over to 
the packing regime, where the isotherm reaches 
a saturation. It is remarkable that the overall 
shape of the isotherm very closely resembles the 
one reported for acetylene [5], including the low- 
pressure feature (see the inset to Figure 1(a) in 
that paper) signalling a very efficient accommo- 
dation. On the other hand, the CO2 molecules 
(which are of similar size as C2H2 and can also be 
viwed as diatomics) do not exhibit this effect. The 
chain length difference could explain the different 
adsorption behavior as a segregation phenomenon 
[16], even without considering a specific adsorp- 
tion. 

Therefore, we have to find out why the chains 
of the same length could behave differently. With 
this purpose we have analyzed the parameters, at 
which the low-pressure condensation occurs. We 
have studied the solution of eq. (10) at to = 2 and 
one-dimensional (ID) (q=2) geometry that mimics 
a channeled host structure, reported in [5]. Inter- 
actions in the neighboring channels were neglected 
and the in-channel interactions were assumed to be 
governed by 32(f)- A fit to the experimental data 
[5] is presented in Figure 3. It seen that the the- 
ory is qualitatively consistent with the data. It is 
found that the low-pressure feature sharpens with 
increasing stickiness A (segment binding energy). 



In order to be free from possible artifacts of the 
mean-field approximation we have studied an ex- 
act solution for ID geometry [17] 
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where J = (3W, H = (3{ijl- W)/2, M = 29-1. It 
is found that there are two regimes (see the inset 
in Figure 3). In a low-activity regime (small A) 
the isotherm exhibits a gradual adsorption with- 
out the low-density condensation, while in the 
high-activity regime (large A) we observe this fea- 
ture. This can explain a striking difference be- 
tween C2H2 and CO2 [5]. Although the latter can 
also be viewed as diatomic (possibly with different 
inter- and intra- molecular constants Q and K) , its 
interaction with the host is much weaker that that 
of C2H2, which is found [5] to form two hydrogen 
bonds with non-coordinated oxygens in the host 
matrix. 

In conclusion, we have found that an interplay 
of the inter- and intramolecular correlations com- 
bined to a variation in the host activity is an ef- 
ficient tool for low-pressure insertion and separa- 
tion of chain molecules. The proposed mechanism 
should be dominant in low-dimensional geometries, 
where the entropic effects (due to the chain confor- 
mation and reorientation) are strongly suppressed. 
We argue that the mechanism is not exclusive to 
particular guest molecules, such that the C2H2 in- 
sertion should be considered as a prototypic ex- 
ample. Therefore, a detailed study of a correla- 
tion between the host geometry, activity and the 
chain length and architecture could open novel per- 
spectives in the storage and separation of other 
technologically important species. In particular, a 
proper control of the guest-induced matrix distor- 
tions [10] or restructuring [11] could lead to a suit- 
able matching or mismatching between the host 
geometry and the guest architecture. This offers 
a promising tool for swithching between the high- 
and low-activity regimes, increasing the insertion 
selectivity in a given domain of controlling ther- 
modynamic parameters (e. g. pressure, density, 
concentration) . 



[1] L. Schlapbach, A. Zuttel, Nature (London), 2001, 
414, 353 

[2] J. S. Slusky et al., Nature (London), 2001, 410, 
343 

[3] O. M. Yaghi et al., Nature (London), 2003, 423, 
705 

[4] G. Ferey et al., Acc. Chem. Res. 2005, 38, 217 



[5] R. Matsuda et al. Nature (London), 2005, 436, 
238 

[6] G. Ferey, Nature (London), 2005, 436, 187 

[7] A. Voronov, I. Luzinov, S. Minko, A. Sidorenko, 

E. Vakarin, M. Holovko, Macromolecules 30, 6929 

(1997) 

[8] A. Voronov, E. Pefferkorn, S. Minko, Macro- 
molecules 31, 6387 (1998) 
[9] S. Minko, A. Voronov, E. Pefferkorn, Langmuir 

2000, 16, 7876 
[10] E. V. Vakarin, J. P. Badiali, M. D. Levi, D. Aur- 

bach, Phys. Rev. B 2001, 63, 014304 
[11] E. V. Vakarin, J. P. Badiali, J. Phys. Chem. B 

2002, 106, 7721 
[12] M. S. Wertheim, J. Chem. Phys. 1987, 87, 7323 
[13] J. P. Badiali, L. Blum, M. L. Rosinberg, 

Chem. Phys. Lett. 129, 149 (1986) 
[14] E. V. Vakarin, Yu. Duda, M. F. Holovko, 

J. Chem. Phys. 107, 5569 (1997) 
[15] J. Chang, S. I. Sandler, J. Chem. Phys. 102, 437 

(1995) 

[16] E. V. Vakarin, J. Chem. Phys. 109, 338 (1998) 



[17] R. J. Baxter, Exactly Solvable Models in Statis- 
tical Mechanics ( Academic Press, NY, 1982) 

FIG. 1. The effective lateral interaction (3W as a 
function of the bulk density r/ at different chain lengths 
m. 

FIG. 2. Mean-field adsorption isotherm for dimers 
m = 2 on a lattice with coordination q = 4. Other 
parameters are A = 800, K = 0.15, Q = 2. The sat- 
uration pressure Po corresponds to 9 — » 1. The inset 
shows typical chain arrangements and correlations in 
different pressure domains. 

FIG. 3. Fitting of the mean-field isotherm (line) to 
the experimental data of Ref.[5] (symbols) for C2H2 
insertion. The saturation pressure Po was taken to be 
5kPa. The inset shows the exact ID isotherms at dif- 
ferent activity regimes (see text). 
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